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FOREWORD 

This  report  contains  the  results  of  an  effort  to  nodlfy  the  loss  reestlaa- 
tlon  procedure  utilized  In  the  UDOSOO  axial  compressor  design  computer  program. 

A  three-sllaenslonal  shock  analysis  was  added  In  order  to  Improve  the  shock  loss 
prediction  accuracy.  The  work  was  done  by  the  Coag>ressor  Research  Group, 
Technology  Branch,  Turbine  Engine  Division,  Aero  Propulsion  Laboratory,  Air 
Force  Wright  Aeronautical  Laboratories  at  Wrlght-Patterson  Air  Force  Base,  Ohio. 
All  debugging  and  code  evaluation  was  done  on  an  IBM  370  system.  The  effort  was 
conducted  by  S.  L.  Puterbaugh  and  Dr  A.  J.  Wennerstrom  from  August  1980  to  August 
1981  under  Project  2307,  Task  SI,  Work  Unit  27,  ''Turbomachlnery  Fluid  Mechanics.” 

The  code  Is  Included  In  the  latest  version  of  DDOlOO  and  Is  available  upon 
request  to  AFWAL/POTX. 
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SECTION  I 


INTRODUCTION 

This  report  describes  the  use  of  a  three-dimensional  shock  model  In  the 
estimation  of  shock  losses  In  a  compressor  rotor  blade  row,  based  on  the  obser¬ 
vation  of  D.  C.  Prince,  Jr.,  In  Reference  1.  In  essence,  the  model  eacplolts 
the  three-dimensional  aerodynamic  sweep  associated  with  a  shock  ^rface  produced 
by  the  spanwlse  distribution  of  a  series  of  normal  shocks  locatM  on  concentric 
stream  surfaces.  The  aerodynamic  sweep  reduces  the  Mach  nuiid>er  normal  to  the 
shock  surface,  hence,  reducing  the  strength  of  the  shock. 

An  Illustrative  example  Is  presented  and  discussed  In  Section  IV  of  the 
report. 
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SECTIOM  II 


BACBSRODHD 

As  stated  In  the  Introduction,  the  shock  aodel  vfalch  was  used  In  the  new 
loss  estlaatlon  procedure  was  discussed  In  a  paper  by  David  C.  Prince,  Jr.  In  the 
Journal  of  Aircraft  (Reference  1).  The  paper  Includes  a  nusiber  of  ezperlaental 
observations  which  Illustrate  two  significant  points  concerning  the  shock  struc¬ 
ture  present  In  txansonlc/supersonlc  blade  rows.  The  two  points  concern  the 
orientation  of  the  passage  shock  which  extends  from  the  leading  edge  of  a  blade 
to  the  suction  surface  of  the  adjacent  blade.  First,  when  the  rotor  is  operating 
at  the  design  point,  the  shock  lies  approximately  axially  (Figure  1);  and  secondly, 
when  the  rotor  Is  operating  at  relatively  high  back  pressures,  the  shock  lies 
approximately  normal  to  the  flow  on  a  stream  surface  (Figure  2).  This  second  case 
represents  operation  at  or  near  peak  efficiency. 

The  consistency  with  which  these  phenomena  occur  allows  for  the  modification 
of  the  existing  shock  loss  calculation  procedure  found  in  UD0300  (Reference  2) 
based  on  their  characteristics.  The  result  Is  a  more  realistic  and,  happily, 
lower  estimation  of  the  loss  Incurred  through  the  presence  of  shocks  in  the  blade 
row.  The  original  method  of  calculation  of  shock  loss  found  in  UD0300  will  now 
be  discussed. 

Losses  are  Incorporated  Into  the  axlsymmetrlc  flow  solution  In  DD0300 
through  relative  total  pressure  loss  coefficients.  The  loss  coefficient  Is  made 
up  of  two  compraents;  the  loss  due  to  diffusion  In  the  profile  boundary  layer, 
and  the  loss  due  to  the  presence  of  shocks. 

The  shock  loss  was  calculated  assuming  a  normal  shock  Is  present  In  the 
blade  passage,  much  like  that  which  exists  when  at  or  near  peak  efficiency.  A 
shock  Mach  number  was  calculated  by  averaging  the  inlet  Mach  number  and  a  suction 
surface  Mach  number  calctilated  by  a  Prandtl-Meyer  expansion  (or  compression) 
through  Input  expansion  (or  compression)  angle  (Figure  2) .  The  total  pressure 
loss  across  the  shock  was  then  calculated  using  standard  normal  shock  parameter 
ratio  equations  at  the  shock  Mach  number. 

This  calculation  produced  a  pessimistic  prediction  at  low  aspect  ratios, 
which  became  even  more  so  at  higher  tip  speeds.  It  \fa8  felt  that  this  was  an 
unrealistic  model  which  severely  limited  attempts  to  explore  higher  tip  speed 
designs. 
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SECTION  III 


THEORY  AND  METHODS 

The  shock  configuration  associated  with  the  naxlauiB  efficiency  operating 
point  was  used  as  the  basis  of  the  model.  This  was  done  due  to  the  geometric 
simplicity  of  that  configuration  In  that  the  shock  Impingement  point  on  the 
suction  surface  could  more  easily  be  determined.  A  comparison  of  the  two  con¬ 
figurations  (Figures  1  and  2)  illustrates  the  rationale  for  this  decision. 

Since  the  efficiency  configuration  was  used,  a  design  which  Is  developed 

utilizing  this  model  must,  obviously,  be  evaluated  In  that  light. 

The  shock  surface  mentioned  previously  Is  created  by  a  spanwlse  distribu¬ 
tion  of  normal  passage  shocks  (Figure  3).  In  low  aspect  ratio  blading,  like 
that  In  Figure  3,  the  aerodynamic  sweep  of  the  shock  surface  at  the  adjacent 
suction  surface  Is  highly  exaggerated  due  to  twist  In  the  blade.  It  Is  the 
combination  of  this  sweep  and  the  leading  edge  sweep  which  Is  used  In  calculat¬ 
ing  the  magnitude  of  the  normal  component  of  the  local  relative  Mach  number 
which  Is  Incorporated  into  the  shock  loss  computation.  This  modified  shock 
Mach  number  reduces  the  total  pressure  ratio  across  the  shock,  thus  reducing 
the  shock  loss  relative  to  the  previous  method. 

The  stream  surface  Intersection  of  two  adjacent  blades  at  each  streamline 
Is  developed  onto  a  flat  surface  normal  to  the  flow  direction  In  order  to  deter¬ 
mine  the  suction  surface  Impingement  point  of  the  shock.  The  impingement  point 
Is  the  Intersection  of  the  suction  surface  of  the  appropriate  blade  and  a  line 
representing  the  shock  which  Is  perpendicular  to  the  inlet  flow  and  intersecting 
the  leading  edge  of  the  adjacent  blade.  The  relative  Mach  number  at  the  shock 
impingement  Is  calculated  by  a  Frandtl-Meyer  expansion  (or  compression)  through  a 
turning  angle  equal  to  the  difference  of  the  inlet  relative  flow  angle  and  the 
suction  surface  angle  at  the  Impingement  point. 

The  aerodynamic  sweep  Is  calculated  for  the  leading  edge  and  the  locus 
of  shock  impingement  points  on  the  suction  surface.  The  sweep  angle,  v.  Is  a 
function  of  several  blade  configuration  parameters  and  Is  given  by: 

V  ■  cos“^  (cos  6  cos  e  sin  (0-Y)-sln  B  sin  e)  -  90* 


where 


B  "  relative  flow  angle 
e  ■  bla<?e  lean  angle 

0  >  streKnllne  slope  angle 

y  ■  computing  station  lean  angle 

The  definition  and  orientation  of  these  angles  are  given  In  Figure  4.  The 
Inlet  and  suction  surface  Mach  numbers  are  then  modified  by  their  associated 
aerod3raamlc  sweep  angles. 

Since  total  pressure  ratio  across  a  normal  shock  Is  not  a  linear  function 
of  Mach  number,  a  simple  three-point  Simpson's  Rule  Integration  of  shock  pres¬ 
sure  ratio  Is  performed  along  the  shock  for  streamlines  at  greater  radii  than 
the  sonic  radius.  Leading  edge,  mid-channel,  and  suction  surface  values  of 
sweep  and  Mach  nuiid>er  are  used  for  the  Integration.  The  mid-channel  value  of 

sweep  and  Mach  number  Is  an  average  of  the  two  end  values.  The  integrated  shock 

pressure  ratio  Is  then  Incorporated  into  the  shock  loss  component  of  the  rela¬ 
tive  total  pressure  loss  coefficient. 

A  finite  value  of  shock  loss  Is  present  In  streamtubes  below  the  sonic 
radius  due  to  a  transonic  bubble  on  the  suction  surface.  The  following  proce¬ 
dure  was  developed  to  address  this  situation.  When  the  sweep  corrected  mid¬ 
channel  Mach  number  Is  greater  than  1.0,  the  shock  loss  Is  calculated  based  on 
the  product  of  the  mid-channel  value  and  the  Inlet  relative  Mach  number.  If 
the  product  Is  less  than  1.0,  the  shock  loss  Is  set  to  0.0. 


SECTION  IV 


OPTION  IMPLEMENTATION 

The  Incorporation  of  the  three-dlnensional  sweep  model  Into  the  loss 
re-estimatlon  calctilatlon  Is  triggered  by  setting  the  value  of  NDEL  at  the 
appropriate  computing  station  less  than  zero.  The  previous  shock  loss  calcu¬ 
lation  procedure  will  be  used  for  positive  values  of  NDEL.  The  only  other 
Input  data  requirement  Is  that  the  number  of  blade  design  passes  must  be 
greater  than  one. 

A  table  of  values  entitled  "Shock  Surface  Sweep  Calculation  Parameters" 
Is  printed  at  the  end  of  the  flow  solution.  The  table  consists  of  streamline 
number.  Inlet  radius,  leading  edge  sweep  angle,  suction  surface  sweep  angle, 
shock  pressure  ratio,  and  calculated  Prandtl-Meyer  expansion  angle. 

The  flow  of  the  program  when  executed  with  the  three-dimensional  shock 
model  Is  as  follows: 

1.  During  the  first  design  pass,  the  losses  are  not  re-estlmated. 

This  Is  done  so  that  a  converged  solution  Is  used  as  a  base  for 
the  re-estlmatlon  passCes). 

2.  The  shock  surface  sweep  values  are  calculated  for  use  In  the 
next  design  pass.  The  calculation  Is  done  In  subroutine  UD0331, 
which  Is  called  from  subroutine  UD03AB. 

3.  After  the  first  design  pass,  the  aerod3nunlc  sweep  angles  calcu¬ 
lated  on  the  previous  pass  are  used  In  determining  the  relative 
total  pressure  loss  coefficients  to  be  used  In  that  pass. 
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SECTICW  V 
TRIAL  EXAMPLES 

In  order  to  substantiate  the  credibility  of  the  geometric  calculations, 
four  test  cases  were  conceived.  Each  was  designed  to  produce  zero  work  and 
each  was  without  any  annular  convergence.  The  Influence  of  blade  lean  and 
leading  edge  meridional  sweep  on  the  aerodynamic  sweep  angle  was  to  be  Inves¬ 
tigated.  The  four  cases  were: 

Case  I  -  zero  LE  sweep  and  zero  LE  lean 

Case  II  -  30"  LE  sweep  and  zero  LE  lean 

Case  III  -  zero  LE  sweep  and  30"  LE  lean  at  the  tip 

Case  IV  -  30"  LE  sweep  and  30"  LE  lean  at  the  tip 

The  results  of  these  example  cases  are  listed  In  Tables  1  thru  4.  All 
trends  of  the  parameter  were  what  was  expected  and,  therefore,  gave  credibility 
to  the  geometric  calculations. 

The  test  of  credibility  of  the  entire  model  was  conducted  In  an  evaluation 
of  the  USAF/AFAPL  HTF  compressor  stage  (Reference  3) .  Since  this  stage  has  been 
fabricated  and  tested,  experimental  results  were  available  which  offered  an  excel¬ 
lent  opportunity  to  evaluate  the  new  shock  loss  model  in  the  light  of  a  design 
result  versus  experimental  result  comparison. 

The  experimental  results  revealed  that  the  Isentroplc  efficiency  curve 
peaked  approximately  5  percentage  points  higher  than  the  design  efficiency. 

When  the  design  was  rerun  with  the  three-dimensional  shock  model  Included,  all 
else  being  equal,  an  Increase  of  approximately  3.75  percentage  points  was  noted. 
Although  this  close  correlation  on  the  basis  of  one  example  could  be  fortuitous. 
It  was  viewed  as  a  highly  promising  result. 
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CONCLUSIONS 

The  shock  loss  calculation  technique  reported  on  herein  appears  to  give  a 
nore  realistic  evaluation  of  the  losses  present  In  a  transonic/supersonic  com¬ 
pressor  rotor  of  low  aspect  ratio.  The  technique  was  evaluated  both  geometri¬ 
cally  and  aerodynamlcally  and  found  to  produce  favorable  results.  The  aerody¬ 
namic  evaluation  proved  most  encouraging,  as  an  analysis  of  the  USAF/AFAPL  RTF 
compressor  predicted  an  efficiency  which  was  1  point  below  experimental  results, 
whereas  the  original  design  utilizing  the  old  shock  model  was  about  5  points  low 
In  efficiency  as  compared  to  the  experimental  results. 

It  must  be  noted  that  because  the  model  Is  based  on  the  shock  configuration 
present  at  the  maximum  efficiency  operating  point,  the  results  must  be  evaluated  In 
that  light. 

The  code  has  been  lAcluded  In  the  most  recent  version  of  UD0300.  The  code 
and  user's  manual  are  available  upon  request  from  AFWAL/POTX. 
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TABLE  1 


RESULTS  OF  TEST  CASE  I 


LEADING  EDGE  VALUES  (DEGREES) 

— 

SL 

BETA 

EPSILON 

PHI 

GAMMA 

NU 

1 

48.02 

0.16 

0.00 

-0.01 

0.12 

2 

49.67 

0.17 

0.10 

-0.01 

0.18 

3 

51.69 

0.18 

0.16 

-0.01 

0.23 

4 

53.88 

0.19 

0.19 

0.00 

0.26 

5 

56.09 

0.19 

0.23 

-0.01 

0.28 

6 

58.22 

0.18 

0.25 

0.00 

0.28 

7 

60.20 

0.18 

0.24 

0.00 

0.27 

8 

62.02 

0.19 

0.21 

0.00 

0.27 

9 

63.67 

0.21 

0.16 

0.00 

0.26 

10 

65.16 

0.27 

0.09 

0.00 

0.28 

11 

66.50 

0.30 

0.00 

0.00 

0.27 

SHOCK  INTERSECTION  VALUES  (DEGREES) 

SL 

BETA 

EPSILON 

PHI 

GAMMA 

NU 

1 

43.80 

-10.98 

0.00 

14.97 

3.18 

2 

45.46 

-  9.40 

0.14 

14.56 

3.60 

3 

47.50 

-  6.91 

0.24 

13.63 

4.28 

4 

49.68 

-  4.57 

0.33 

12.64 

4.90 

5 

51.90 

-  2.41 

0.38 

11.49 

5.42 

6 

54.02 

-  0.56 

0.41 

10.33 

5.84 

7 

56.01 

0.99 

0.40 

9.22 

6.17 

8 

57.84 

2.28 

0.35 

8.20 

6.45 

9 

59.48 

3.35 

0.27 

7.32 

6.71 

10 

60.94 

4.25 

0.15 

6.47 

6.91 

11 

62.22 

4.97 

0.00 

6.04 

7.19 

SL  ■  Streamline  Humber 

BETA  *  Relative  Flow  Angle 

EPSILON  ■  Blade  Lean  Angle 

PHI  >  Stream  Surface  Pitch  Angle 

GAMMA  *  Lean'  Angle  In  the  Meridional  Plane 

HU  -  Aerodynamic  Sweep  Angle 
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TABLE  2 


RESULTS  OF  TEST  CASE  II 


LEADING  EDGE  VALUES  (DEGREES) 

SL 

BETA 

EPSILON 

PHI 

GAMMA 

NU 

1 

48.51 

2.82 

0.00 

30.00 

21.26 

2 

50.08 

2.76 

0.27 

30.00 

20.80 

3 

52.00 

2.74 

0.49 

29.97 

20.14 

4 

54.00 

2.64 

0.64 

30.10 

19.36 

5 

56.18 

2.57 

0.70 

29.95 

18.40 

6 

58.19 

2.52 

0.69 

29.91 

17.47 

7 

60.05 

2.47 

0.63 

30.08 

16.67 

8 

61.73 

2.40 

0.51 

29.98 

15.77 

9 

63.21 

2.29 

0.36 

30.01 

14.98 

10 

64.50 

2.15 

0.19  , 

30.00 

14.22 

11 

65.61 

2.00 

0.00 

30.00 

13.53 

SHOCK  INTERSECTION  VALUES  (DEGREES) 

SL 

BETA 

EPSIUEf 

PHI 

GAMMA 

NU 

1 

44.96 

-13.33 

0.00 

39.73 

18.27 

2 

46.49 

-11.25 

0.20 

39.60 

18.82 

3 

48.37 

-  7.98 

0.39 

39.21 

19.82 

4 

50.35 

-  4.92 

0.56 

38.82 

20.62 

5 

52.35 

-  2.08 

0.68 

38.12 

21.11 

6 

54.25 

0.04 

0.69 

37.50 

21.42 

7 

56.03 

2.34 

0.63 

37.05 

21.60 

8 

57.66 

3.87 

0.51 

36.42 

21.49 

9 

59.12 

5.01 

0.36 

35.96 

21.30 

10 

60.44 

5.80 

0.19 

35.49 

20.96 

11 

61.63 

6.49 

0.00 

35.25 

20.72 
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TABLE  3 


RESULTS  07  TEST  CASE  III 


LEADING  EDGE  VALUES  (DEGREES) 

SL 

BETA 

EPSILCm 

PHI 

GAMMA 

1 

48.02 

30.53 

0.00 

-0.01 

22.19 

2 

50.11 

28.71 

0.10 

-0.01 

21.69 

3 

52.59 

26.90 

0.16 

-0.01 

21.15 

4 

55.18 

25.13 

0.20 

-0.00 

20.51 

5 

57.71 

24.45 

0.23 

-0.01 

19.78 

6 

60.06 

21.91 

0.25 

0.00 

18.98 

7 

62.20 

20.56 

0.24 

0.00 

18.21 

8 

64.11 

19.29 

0.21 

0.00 

17.38 

9 

65.81 

18.14 

0.16 

0.00 

16.57 

10 

67.31 

17.35 

0.09 

0.00 

16.00 

11 

68.64 

16.54 

0.00 

0.00 

15.37 

SHOCK  INTERSECTION  VALUES  (DEGREES) 

SL 

BETA 

EPSILON 

PHI 

1 

43.80 

18.65 

0.00 

19.48 

25.97 

2 

45.64 

22.80 

0.13 

18.58 

28.13 

3 

47.78 

28.27 

0.24 

16.61 

30.70 

4 

49.98 

33.38 

0.32 

14.72 

33.40 

5 

52.10 

38.05 

0.38 

12.79 

36.05 

6 

54.09 

42.08 

0.41 

11.03 

38.54 

7 

55.91 

45.54 

0.40 

9.48 

40.86 

8 

57.54 

48.40 

0.35 

8.17 

42.90 

9 

59.01 

50.76 

0.27 

7.10 

44.71 

10 

60.30 

52.96 

0.15 

6.16 

46.44 

11 

61.42 

54.69 

0.00 

5.68 

47.97 

1< 


TABLE  4 


RESULTS  OF  TEST  CASE  IV 


leading  edge  values  (DEGREES) 

SL 

BETA 

EPSILON 

PHI 

GAMMA 

NU 

1 

48.51 

35.71 

0.00 

30.00 

41.67 

2 

50.54 

33.86 

0.27 

30.00 

40.78 

3 

52.93 

32.05 

0.49 

29.97 

39.64 

4 

55.42 

30.25 

0.63 

30.09 

38.28 

5 

57.84 

28.49 

0.70 

29.95 

36.63 

6 

60.07 

26.89 

0.69 

29.91 

34.97 

7 

62.07 

25.52 

0.63 

30.08 

33.46 

8 

63.82 

24.13 

0.51 

29.98 

31.83 

9 

65.34 

22.84 

0.35 

30.01 

30.32 

10 

66.63 

21.84 

0.18 

30.00 

29.05 

11 

67.73 

20.79 

0.00 

30.00 

27.79 

SHOCK  INTERSECTION  VALUES  (DEGREES) 

SL 

BETA 

EPSILON 

PHI 

GAMMA 

NU 

1 

44.96 

22.61 

0.00 

42.47 

40.43 

2 

46.68 

27.83 

0.20 

42.05 

42.28 

3 

48.66 

34.38 

0.39 

41.04 

44.70 

4 

50.65 

40.16 

0.56 

40.12 

47.11 

5 

52.56 

45.11 

0.68 

38.98 

49.40 

6 

54.31 

49.18 

0.69 

37.99 

51.53 

7 

55.91 

52.58 

0.63 

37.28 

53.49 

8 

57.34 

55.19 

0.51 

36.47 

55.12 

9 

58.64 

57.24 

0.35 

35.88 

56.51 

10 

59.80 

58.98 

0.19 

35.32 

57.75 

11 

60.84 

60.39 

0.00 

35.06 

58.83 
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appendix  a 

DETAILED  DERIVATION  OF  AERODTNAMIC  SWEEP  ANGLE 

V  defined  as  the  angle  between  Fa  and  B  alnus  90* 

¥  -  cos  e  sin  (0-y)T  -  sin  e  J  +  cos  e  cos  (0-Y)k 
Fa  ■  cos  BT  +  sin  $j 

^s  *  B  ■  |Fg||B|cos  (90fr) 

But  |Fa|-|B|  -  1 

*.  V  -  cos”*(|Fa||B|)  -  90 

Fa*B  -  (cos  B1  +  sin  bT)*(cos  €  sln(0-Y)l  -  sin  ej  +  cos  e  cos  (I-yW 
)Fa*B|«  cos  B  cos  e  sin  (0-Y)  -  sin  B  sin  e 

V  -  cos”*  (cos  B  cos  e  sin  (0-y)  -  sin  B  sin  e)  -  90 
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